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Abstract. We show that, for aU odd natural numbers A'', the A'-torsion 
points on an elhptic curve may be placed in an TV x TV grid such that 
the sum of each column and each row is the point at infinity. 



1. Introduction 

Let be a positive integer, and consider the integers 1, 2, . . . , N"^. An 
N X N grid containing these consecutive integers such that the sum of each 
column and each row is the same is cahed a magic square. (This is usuahy 
cahed a semi-magic square in the hterature; see [6J.) For example, when 

= 3, we have the grid 
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where the sum of each column and each row is 15. 

We need not limit ourselves to a grid with integer entries. The author 
of [l], inspired by the discussion in \T, §1.4], considered the problem of 
arranging the 9 points of inflection on an elliptic curve in a 3 x 3 magic 
square. That is, it is possible to arrange the points of order 3 in a 3 x 3 grid 
so that the sum of each row and each column is the same, namely the point 
at infinity. We generalize this result: 

Theorem 1. Let N > 1 be an odd integer, let E he an elliptic curve defined 
over an algebraically closed field with characteristic not dividing N . Then 
the N"^ points of order N on E can he placed in an N x N magic square 
such that the sum of each column and each row is the point at infinity O. 

We construct such a grid using Lehmer's Uniform Step Method, as motivated 
by the discussion in [4]. In particular, the theorem holds for any group G 
such that the A^-torsion G[N] ~ (Z/A^Z) x (Z/NZ). 

The author would like to thank Jiu-Kang Yu for helpful comments. 

2. Semi-Magic Squares over Abelian Groups 

As stated above, we define a magic square to be an A x A grid containing 
the consecutive integers 1 through A^ such that the sum of each column and 
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each row is the same. Strictly speaking, this is a semi-magic square, but we 
abuse notation shghtly for the sake of brevity. We do not hmit ourselves to 
constructing magic squares with integer entries. Indeed, we will construct 
an N X N magic square for a certain class of abelian groups. 

Let G be an abelian group under ©. Given P £ G, denote [— 1]^^ as its 
inverse and [0]P = O as the identity. For each nonzero integer m, denote 
[m]P as [±1]^^ added to itself \m\ times, where "it" is chosen as the sign 
of m. Denote G[m] C G as that subgroup consisting of points P G G such 
that [m]P = O. We will always assume that G is chosen such that for some 
positive integer there is a group isomorphism 

(2.1) V : (Z/iVZ) X (Z/NZ) — ^ G[N]. 

We have a bijection {1, 2, . . . , A^^} ^ (Z/NZ) x {Z/NZ) given by 

(2.2) (p: (^k-l (mod iV), 

where [-J is the greatest integer function. That is, if 1 < A; < A^^ then we 
can write k — 1 = m + N n for some unique < m, n < N , and so we map 
k I— > (m, n). This means we have a bijection of sets 

^Pocj): {l,2,...,iV2} G[N]. 

We will use this identification to place the elements in G[N] in an N x N 
magic square. 

There are two examples in particular which will be of interest to us. Upon 
fixing A^, the group G = (Z/NZ) x (Z/NZ) satisfies the criterion above. As 
another example, fix an algebraically closed field F and let E be an elliptic 
curve defined over F. We may choose G = E(F) as the F-rational points on 
E, where we have a non-canonical isomorphism G[N] ~ (Z/NZ) x (Z/NZ) 
only when the characteristic of F does not divide A'^. (For more properties 
of elliptic curves, see [3].) 

3. Uniform Step Method 
Fix a positive integer A^. Let G be an abelian group under ©, and assume 
G[N] = {Ri, R2, Rk, Rn'2} - (Z/NZ) X (Z/NZ) . 

We wish to place these A^^ elements in an A^ x A^ grid such that the sum of 
each row and the sum of each column, as an element in G, is the same. We 
use an idea of D. H. Lehmer from 1929, known as the Uniform Step Method. 
To this end, we are motivated by the discussion in [4, Chapter 4]. 

Given an N x N grid, we consider its entries in cartesian coordinates. For 
the moment, fix integers a, b, c, and d, and consider placing the element 
Rk € G'[A^] in the (xk,yk) position. After arbitrarily placing Ri in the 



k-1 
N 



(mod n] 
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(xi, yi)-position, we will define Xk and yk by the recursive sequence 

k-1 



Xk = xi + a {k — 1) + b 
yk=yi+c{k-l) + d 



N 
k-1 

N 



(mod A^) 

for l<k< N"^. 

(mod N) 



We will exhibit conditions on these integers a, b, c, and d such that the 
sequences above indeed generate a magic square. 

Proposition 2. If N is odd and relatively prime to {ad — be), then the 
sequence {x^, yk) places exactly one Rk in each of the N'^ cells of the N x N 
grid. 

Proof. It suffices to show that (xfcj,yfcj = {xk^^y^^) only when /ci = /c2; for 
then we would have N"^ different points so they must fill in the entire grid. 
Using the bijection (p as in ()2.2p note that we may write 

Xk = xi + am + bn (mod N) 

where (m,n) = 0(A;). 

Vk = yi + cm + dn (mod rs ) 
Write (mi,ni) = 4'{ki) and (7712,^2) = 4>{k2), so that 

a (mi — 171,2) + b (ni — 712) = (mod N) 
( ki,yki) ( ki,yk2) g _ ^^2) + d (r7i — 712) = (mod A^) 

Since ad — be (mod N) is invertible, we see that this happens if and only if 

(j){ki) = {mi, ni) = (7772, ^2) = 0(A;2) 
and so fci = /c2- D 

Proposition 3. // is relatively prime to a and b, then the sum of the 

entries in the ith column is O. If N is relatively prime to c and d, then the 
sum of the entries in the jth row is O. 

Proof. The entries in the ith column consist of those corresponding to k 
such that Xk = i. Similarly, the entries in the jth row consist of those Rk 
corresponding to k such that yk = j- Hence the sum of the entries in the 
ith column and jth row are 



Rk and Rk, respectively. 



a;fe=« yk=j 



First we determine the values of k which occur in the ith column. Since N 
is relatively prime to a and b, there are exactly N pairs {m,n) G {Z/NZ) x 
(Z/A^Z) satisfying am + bn = i — xi (mod A^); indeed, given any 777 we 
can solve for 77, and vice-versa. Hence there are exactly A'^ integers k = 
1 + 777 + A^ 77 (mod A^^) such that Xk = i, which we denote by ka- If we 
denote {ma,na) = ^(^o) using the bijection in (12. 2p . then it is clear we have 

{. . . , 777a, ...} = {..., 77a, ...} = 'L/N'L. 
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Now we compute the sum of the values m the ith column. Using the 
group isomorphism in (|2.ip . denote P = ip {{1,0)) and Q = Tp {{0,1)) so that 
we have Rk = [m]P [n]Q when {m,n) = <j){k). This gives the sum 



^ -Rfe - ^ Rka 
Xf^=i a 



E 



{[ma]P e K]Q) = [m']P e [n']Q 



where we have set 

E 



m 



n 



rrir, 



m 



N{N -1) 



(mod N). 



Since is assumed odd, this sum is a multiple of N so that [m']P = [n']Q = 
O. Hence the sum of the entries in the iih. column is indeed O. 

A similar argument works for the jth row. □ 

We summarize this as follows. 

Theorem 4. Let G he an abelian group under ©, and assume that there is 
a positive odd integer N such that 

G[N] = {Ri, R2, Rk, Rm} ^ {Z/NZ) X {Z/NZ) . 

Fix integers a, b, c, and d relatively prime to N such that {ad — be) is also 
relatively prime to N, and consider the sequence {xk,yk) defined by 



Xk = xi + a {k — 1) + b 
yk = yi+c{k-l) + d 



k-1 

N 
k - 1 

N 



(mod N) 
(mod N) 



for l<k< N"^. 



The N X N grid formed by placing Rk in the {xk,yk) position is a magic 
square, where the sum of each column and each row is the identity O. 

We remark that this method does not exhaust all ways in which a magic 
square can be generated. For example, this method does not seem to work 
for even. Indeed, the sum of each column and each row involves the 
expression A^(A^ — l)/2, which in general is not a multiple of N . Also, when 

= 4, we have the magic square 
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It is easy to check that such a square cannot be generated by a sequence 
{xk,yk) for any o, b, c, or d. This first appeared in 1514 in an engraving by 
Albrecht Diirer entitled "Melencolia." 
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4. Applications 

We can specialize a, 6, c, and d to generate examples of magic squares. 

Corollary 5. Let G be an abelian group under ©, and assume that there is 
an odd positive integer N such that 

G[N] = {Ri, R2, . . . , Rk, . . . , Rn^} ^ (Z/iVZ) X (Z/NZ) . 

Then these elements can be placed in an N x N magic square such that the 
sum of each column and each row is the identity O. 

Proof. We follow the construction using a method first outlined by De la 
Loubere in 1693. (An example of how this method works follows at the end 
of the paper.) Using Theorem [H set a = 1, 6 = c = —1, and d = 2. As 
is odd, it is relatively prime to these integers as well as the determinant 
ad — be = 1. □ 

Remark. Theorem [1] follows from this corollary, since the group E[N] of 
A'^-torsion points on an elliptic curve E is isomorphic to (Z/A^Z) x (Z/A^Z). 

The following was pointed out to the author by J.-K. Yu. Upon choosing 
the basis {P, Q} for G[N] given by P = V'((1,0)) and Q = ip{iO,l)), we 
may write Rk = [m]P © [n]Q when {m,n) = (p{k). (Here, we use the maps 
defined in (j2.ip and (j2.2p .) In this way, we may identify R^ with {m,n). If 
we choose a = d = 1 and 6 = c = 0, then we have a magic square upon 
placing {m,n) = (j){k) in the (x^, yfc)-position. In general, if for odd A^ we 
have an N X N Latin Square with the (m, n)-position having entry amn then 
we may place (m, amn) in the (x^, yfc)-position. (For more on Latin squares, 
see [5].) 

We discuss a specific example by considering the 3-torsion on elliptic 
curves; to this end, set A^ = 3. We explain how this construction generalizes 
that in [1]. Consider an elliptic curve defined over the complex numbers C, 
and let G = E{C) be the group of complex points on the curve. Then it is 
well-known that we can express the 3-torsion as 

E[3] = {A, B, C, D, [-l]A, [-1]B, [-1]C, [-1]D, O} ~ (Z/3Z) x (Z/3Z) 
where B = A (B D and [—1]B = G(BD. If we label these points as 



Ri 


= 0, 


R4 


= D, 


R7 = [-l]D 


R2 


= [-l]B, 


R5 




Rs = C, 


R3 


= B, 


Rq 




Rg = A; 



then we can use the magic square from the introduction to place the 3-torsion 
in a magic square: 
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We can also compute this magic square using the method in the proof of 
the corollary. Choosing the basis P = [— 1]-B and Q = D; \i can be easily 
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checked that Rk = [m]P © [n]Q when (m,n) = (pik). If we also choose 
(xi,yi) = (2,2) as the center of the 3x3 grid, then Rk may be placed in 
the (xfc, y/c)-position, where 

k — 1 

(mod N) 

for 1 < A; < N"^. 

(mod A^) 



Xk = xi + {k- 1) 



yk = yi-{k-l) + 2 



N 
k- 



1 



N 



As mentioned before, this is known as De la Loubere's method or the Siamese 
method. 
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